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PRODUCTS 

S. CAROLINA GARCIA-MARTINEZ AND JONATAN HERRERA 


Abstract. In this paper we study the local rigidity of metrics defined on a compact manifold 
M with boundary dM and satisfying both constant scalar curvature on M and constant mean 
curvature in dM. We present some geometrical hypotheses ensuring local rigidity first for the 
general Riemannian case, and then, for the warped metric one. These conditions arise from 
the study of a spectral problem which is not included within the classical problems (Neumann, 
Steklov,...) and that we call “mixed eigenvalue problem”. Finally, we apply our results for the 
spatial slice of the Anti-de Sitter spacetimes. 
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1. Introduction 

The classical Yamabe problem consists in showing that every Riemannian compact manifold, 
without boundary, admits a conformally related metric with constant scalar curvature. In 
dimension two, such a problem follows from the Uniformization Theorem of Riemann Surfaces. 
For higher dimensions, the problem was formulated by Yamabe in He], where the first steps 
towards a proof of the existence of a solution were given. After the combined efforts of Trudinger 
[24] . Aubin [5] and Schoen [22], the Yamabe problem was completely solved. It is important 
to highlight that constant scalar curvature metrics can be characterized variationally as critical 
points of the Hilbert-Einstein functional in conformal classes. It is known that the minimum 
of this functional in a conformal class is unique, see [2]. However, in many cases a rich variety 
of constant scalar curvature metrics arise as critical points that are not necessarily minimizers, 
and for this reason, it is very interesting to find conformal classes where the Yamabe problem 
has multiple solutions. Among others, bifurcation techniques can be applied in order to obtain 
multiplicity results. The existence of multiple solutions for the Yamabe problem in different 
settings has long been studied in the literature, see for instance HUSUZUE]. 

Now, if the compact manifold considered in the Yamabe Problem has a nonempty boundary, 
several possible boundary conditions can be studied. For instance, from the point of view of 
conformal geometry, a geometrical condition would have to involve the mean curvature of the 
boundary. This important observation led Escobar in 1992 (see m) to study the problem of 
finding smooth metrics with constant scalar curvature and minimal boundary inside a given 
conformal class. In a natural continuation of the latter, the same author in m analyzed 
the conformal deformation of a metric to a scalar flat metric with constant mean curvature 
on the boundary generalizing the famous Riemann Mapping Theorem to higher dimensions. 

l 
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Subsequently, he also obtained results for the Yamabe problem with boundary under mixed 
constraints (see my- 

(1) aV + bA = l 

where V and A are the volume of manifold and the area of its boundary, respectively, and with 
a > 0 and b real numbers. More precisely, he studied a problem of existence of conformal metrics 
with constant scalar curvature R{g ) in the manifold and constant mean curvature H(g) on the 
boundary, which are related by 


(n — 1 )bR(g) = 2 nH(g)a. 


The solutions of this problem are critical points of the functional formulated by Gibbons-Hawking 
and York (13 , 20] (which we will refer to as the GHY-functional ) restricted to the constraint 0 - 
The existence of such critical points was proved by Escobar for manifolds of nonpositive type 
(i.e., the first eigenvalue of the conformal Laplacian with zero boundary condition is nonpositive) 
and for almost any manifold of positive type if b is sufficiently small. In m and na, Han and 
Li also showed the existence of solutions when a > 0, the manifold is of positive type and either 
is locally conformally flat with umbilic boundary, or n > 5 and the boundary has a non-umbilic 
point. 


Taking advantage of the variational approach, our aim in this paper is to study the local 
rigidity of families of metrics defined on a fixed compact manifold M with boundary dM and 
satisfying both constant scalar curvature in M and constant mean curvature on dM. Roughly 
speaking, a family T of metrics is (locally) rigid if given another metric g solution of the Yamabe 
problem and sufficiently close to some element of J-, then g belongs to T . 

For our results, we will give conditions to ensure the non-singularity of the second variation 
of the GHY-functional under the constraint ([!]), which will ensure the rigidity as a consequence 
of the Implicit Function Theorem (see m Appendix A]). Such a conditions will be obtained by 
making lower estimates of the first eigenvalue of a particular kind of spectral problem, which we 
have called mixed eigenvalue problem (see Section [3] for a general background on this problem). 

When we consider warped product spaces with constant scalar curvature, a special family 
of metrics appears. As it is well known, warped spaces present a natural foliation given by 
hypersurfaces with constant mean curvature (its fibres). Under the assumption that certain 
initial fibre of the space is minimal, the foliation is naturally identified with an one-parametric 
family of metrics satisfying previous conditions, and so, susceptible for our results on rigidity. 
Such results can be interpreted as local uniqueness of the foliation from a metric viewpoint. 

As an application of the above, we prove that three-dimensional warped products with rnono- 
tonically increasing warping function are always globally rigid (see Theorem 5.3). This result is 
specially interesting when the fibre of the warped product is a sphere since it becomes the spatial 
slice of the Anti-de Sitter Schwarzschild spacetime, a classical relativistic model of a spacetime 
containing a black hole. Our result allows us to deduce the local uniqueness of the model under 
metric variations preserving both the causal structure (i.e., its conformal structure) and the 
boundary of the black hole (see [25] for a background on general relativity and the Anti-de 
Sitter model). 


This work is organized as follows. In Section [2j besides to fixing the notation, we include 
some preliminaries on bifurcation theory, as well as the criterion that we will use to prove 
local rigidity. In Section [3j we discuss a general background about the spectrum of the mixed 
eigenvalue problem, including classical results as the Courant’s nodal Theorem and the Rayleigh 
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characterization for eigenvalues (which will be necessary in the sequel). Moreover, we present 
some results for the particular case of warped spaces. 

In Section [4j we first introduce the variational setting for the Yamabe problem with boundary, 
as we described before. Then, we obtain our main results about the stability of the solutions of 
the Yamabe problem for the general Riemannian case (Theorem 4.8). The case of warped spaces 
is studied in Section [5] where first we introduce the one-parametric family of metrics associated 
to the foliation ( |18[ ), and then we obtain the main result for the warped case (Corollary |5.2[ ). Ad¬ 
ditionally, in Section 5.1 we present several examples where our results are applicable, including 
the case of the spatial fibre of the Anti-de Sitter Schwarzschild model. We have also included 
in Section |5.2| a general procedure to study both the metric rigidity and bifurcation in warped 
spaces not covered by our general results. Finally, in the Appendix (Section [6]), we compute 
the first and second variation of the GHY-functional. This computation has been performed by 
several authors (see [3| for instance), but we include it here for the sake of completeness. 


2. Preliminaries 

The purpose of this section is to state the basic elements, results and notation that we are going 
to use in the rest of the paper. Let ( M n ,g ) be an arbitrary n-dimensional Riemannian compact 
manifold with non-empty boundary dM and n > 3. We will assume that the boundary dM 
is an (n — l)-dimensional smooth manifold (i.e., a hypersurface) formed by different connected 
components. Additionally, we will also assume that such a components are grouped in two 
disjoint sets dM = Si U S 2 (see for instance Figure [lj. As we will see in subsequent sections, 
our intention here is to group all the components sharing the same (constant) mean curvature, 
allowing us to work with them as a single component. 



For a given system of local coordinates x\, ■ ■ ■ ,x n around a point p, the metric g and the 
volume element will be represented as 

g = gij dx l dx\ dv(g) = \J det (gtj) dx 1 A • • • A dx n , 
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with indexes varying from 1 < i,j < n. When p belongs to dM, we will assume that the 
coordinate x n is normal to the hypersurface dM and pointing inward. In particular, the area 
element restricted to the boundary will take the following form: 

da(g) = yjdet(g a p) dx 1 A • • • A dx n_1 

where 1 < ck,/3 < n — 1. [j] Finally, the scalar curvature in ( M,g ) will be denoted by R(g), while 
the mean curvature on each H*. will be denoted by H{gif) for k = 1, 2. 


2.1. Bifurcation theory on metric variations. Here we introduce some basic framework 
about bifurcation theory that we will need along the rest of the paper. We refer the reader to 
HU [23] and the references therein for general background. 

As a first step to define the concepts of rigidity and bifurcation, we have to define a natural 
norm on the space of Riemannian metrics. For this, let us denote by S k (M ), with k > 2, the 
space of all symmetric (0, 2)-tensors of class C k defined on M. The space of Riemannian metrics 
on M, denoted here by A4, is a subspace of S k (M) with the structure of an open cone and 
tangent space T g M naturally identifiable with the entire S k (M). Now, consider an auxiliary 
Riemannian metric gA in M defining both a connection Vyi and an inner product {{■,■)) a on 
the space S k (M). Then, we define the norm || • || C jt given by 


\\T\\c k = maxj = i v .. ifc max peM ||V^ ) T(p )|| J 4 
For g € M, let us denote 


(2) Conf(g) = {fg : f G l~L l {M) and / > 0} 

the space of 'H 1 (M)-conformal metrics of g , which is clearly identified as an open subset of 
l~L l {M). We will consider in Conf(g) the differential structure induced by R l {M). 

Now, let us consider an one-parameter family of metrics {<7 a}a£/ C M, where I is an arbitrary 
open interval on M. We will assume that, for all A, the metric g\ has: 

(a) constant scalar curvature R(g\) in M, 

(b) constantly zero mean curvature on Hi (i.e., Hi is a minimal hypersurface) and 

(c) constant mean curvature (CMC for short) H(g 2 ) on H 2 . 

In this setting, A* is said to be a point of bifurcation if there exist a sequence {A n } n C I and a 
sequence {g n }n C M. satisfying: 

(1) lim A n = A* and lim g n = g\ t (the latter with the topology induced by || • ||( 7 fc), 

n —¥ 00 n —¥ 00 

(2) For all n, g n determines constant scalar curvature in M; and makes Hi a minimal 
hypersurface and H 2 a CMC-hypersurface. Moreover, the scalar curvature in M and the 
constant mean curvature on H 2 are R(g\ n ) and H((g\ n ) 2 ), respectively, 

(3) For all n, g n belongs to the R 1 (Af)-conformal class of g\ n , but g n / g\ n . 


^As a convention, we will always assume that, when we work in coordinates, the indexes i,j will vary between 
1, • • • , n, while the indexes a , /? will vary in 1, • • • , n — 1. 

2 Here, ?^ 1 (M) is the Sobolev space of L 2 functions on M with first derivatives in L 2 . 
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If A* is not a bifurcation point, we will say that the family {<?a}a is locally rigid at A*. 

In the literature, we can find several criteria for the existence of bifurcation points as well as 
for (local) rigidity. In general, such results require a variational problem whose critical points 
are the objects of interest. For our particular case, given a family {#a}a of metrics, we need a 
path of C fc -functionals (with k > 2) F x : Conf(^A) —> M whose critical points are metrics g in the 
conformal class of g\ with constant scalar curvature R(g) = R(g\) for a fixed A E J; and making 
Si minimal and S 2 of constant mean curvature with H(g 2 ) = H({g\) 2 )- Then, if the metric g\ t 
is a non-degenerate critical point for some A*, the fiber bundle version of the Implicit Function 
Theorem ensures the local rigidity on A* (see [ 161 Appendix] for details). For bifurcation, in 
addition to the degeneracy of g \ t , we also require a variation in the associated Morse index to 
the critical points in a close neighbourhood of A*. Namely, let us denote by i(F,g) the Morse 
index of g, that is, the dimension of the maximal subspace of the tangent space T g .M. where the 
second variation 5 2 {F) g is negative definite. Then, A* is a bifurcation point if for any Ai, A 2 close 
enough to A* with Ai < A* < A 2 , i(F Xl ,g aJ / i(F X2 ,g\ 2 ) (see 02E Theorem A. 2] for details). 

In this work, we will focus our attention on local rigidity, even so our studies allow us to 
do some analysis involving bifurcation points (see Section [5.2[ ). In particular, the desired non¬ 
degeneracy will follow by showing that all the eigenvalues of the associated Jacobi operator 
(which is diagonalizable under our conditions) are positive. In this case, the spectrum of the 
Jacobi operator lead us to a mixed eigenvalue problem. 

3. Spectrum of the mixed eigenvalue problem 

Our aim in this section is to give a general background of the so-called mixed eigenvalue 
problem, giving the basic properties that we will need through this paper. This kind of problem 
cannot be consider a Neumann problem nor a Steklov one but a mix of them, can be described 
in the following way 


A g f = (G + /3) f in M, 

- 7 ^/ = {J + J) f on dM 

where A g is the Laplace-Beltrami operator with nonnegative spectruir0 G : Ai —> M and 
J : dM —» M are functions; d/dn denotes the inward normal derivative; and / 6 Jf 1 (M). 

To what extent are the classical procedures for the study of, say, Dirichlet problems applicable 
(or, at least, adaptable) to the mixed one? As a first step to answer this question, it is essential 
to describe the mixed problems variationally. Let us define two bilinear forms V and £ in the 
following way: 


( 3 ) 


T ) (<£>,V’) = / [g(grad(<p),grad(il)))-G(pilj]dv(g)- JLp^da^g) 


im 


/ dM 


/ <pil>dv(g)+ / (pip da (g). 

J M J dM 

Then, for a fixed /3 G M, we define the functional 


3 ; 


i.e., A g = — divg(grad). 
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= T>i}p) - fi £ fip) 

where T>((p) := T>(ip,cp) and £(<p) := £(ip,ip)] and observe that the first variation of such a 
functional becomes: 

WO = 2 _ P £ VO) • 

Now, the first Green identity leads us to: 


£>(¥>> VO - /?£(<£, VO = [ fi {A g p - Gip - Pip) dv(g) + [ p (-d n ip - Jip - Pip) dcr(g). 

Hence, ip is a critical point for the functional J 7 , if and only if it is a solution for the mixed 
eigenvalue problem ([3]). By using this variational approach, we can re-obtain some classical 
and well-known results for the mixed case. For instance, the Courant’s Nodal Theorem follows 
directly (see [8] Page 452]) as well as the classical characterization for the eigenvalues was 
obtained by Courant [8] and Rayleigh [21] . In fact, it follows that the eigenvalues are determined 
by a sequence {/3j}ieN) repeated according to their multiplicity, and such that lim^oo p i = oo. 
Moreover, 


fin = 


mm 


2%) 


£{<P) 

where each ipi (with 1 < i < n — 1) is the eigenfunction associated to P i and 

{ipi,..., Vn-i} L = w G H l (M) : £(ip, (Pi) = 0, Vi = 1,..., n — 1}. 
For convenience, when n = 1 the minimum is taken on the whole 
Now, observe that the first Green identity lead us to: 




\grad(ip)\ 2 g dv(g) = 


IM 


ipAgipdv(g) - / ipd n ipda(g) 


IdM 


j (p + G) <p 2 dv(g) + [ (fi + J) <p 2 da(g) 
J M JdM 


= fi[ V dv (9) + / V dcr(g) 
\Jm JdM 

+ [ Gip 2 dv(g) + f Jtp 2 da(g). 
J M JdM 

Then, if we assume the following normalization 


( 4 ) 


/ l P 2 dv(g) + / ip 2 da(g) = £(ip,ip) = 1 , 
Im JdM 


any eigenvalue p should satisfy 


(5) P= \grad(ip)\ 2 dv{g) — / Gip 2 dv(g) - / Jip 2 da(g). 

Jm Jm JdM 

In particular, we can obtain easily the following simple lower estimate for the eigenvalues 
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Proposition 3.1. The eigenvalues (3 satisfies 

P>-{G + + J+) . 

where G + = maxjvf{0, G} and J + = maxg^iO, </}. 

Proof: First, let us recall that G (and analogously J) can be expressed as G = G + — G~ 
where G + = maxM{0,G} and G~ = maxjv/{0, —G}. Thus, from ([5]) we have that 

P = [ \grad(ip)\ 2 dv(g) - [ Gip 2 dv(g) - [ Jp 2 da{g)> 

JM JM JdM 

> —G + f ip 2 dv(g) — J + f p 2 da(g). 

JM JdM 

Using now the normalization in Q, we have that both previous integrals are less than or equal 
to 1, so we get 

P > - (G+ + J + ) . 

□ 

Previous lower bound of the eigenvalues will be clearly insufficient to obtain rigidity results 
in the forthcoming sections. However, for numerical approaches, it determines a start-point for 
algorithms looking for the first eigenvalue. 

3.1. The mixed eigenvalue problem for warped metric spaces. In this section we will 
focus our attention on the warped product spaces, that is, an n -dimensional Riemannian space 
( M,g ) where 

M = (n, rfi) x P, g = dr 2 + a 2 (r)g p , 

(P, g p ) is an (n — l)-dimensional closed Riemannian manifold and a is a C k positive function 
with k > 2 . Along this section, we will also assume that the previous functions G and J 
only depends on the parameter r. With these assumptions, the mixed eigenvalue problem can 
be reduced to a Sturm-Liouville-type problem with boundary conditions. In fact, when g is a 
warped metric, the Laplace-Beltrami operator splits as 

Agip = + “2 Ap 

where A p denotes the Laplace-Beltrami operator on ( P,g p ). Therefore, the first equation in 
© becomes 

(6) - jd r (a^drip) + ^.A P (p = (G + J) 

Next, by separation of variables, we assume that a solution <p : M —> M of the above problem 
can be split as, 

<p(r,x) = ipm(r)(p P (x) 

where (p^_ : (ri,^) —> M and tpp : P —> M. Moreover, if we assume that tpp = p l p is a non-zero 
eigenfunction for A p associated to an eigenvalue fii then ([6]) can be written as follows 

~~Jrr: dr {u n ~ l drm) + = (g + p) m 
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which is a Sturm-Liouville equation on the interval (ri,r 2 ). Hence, joining previous equation 
with the initial conditions on ( 3]) we obtain, for each eigenvalue Pi of A p, the following problem^ 


(a n 1 drm) + ^2 


(7) 


= (G + p) ip R , 


-<Pw(ri)= (J(n) + 0) <pm(ri), 
m(r 2 ) = (J(r 2 ) + P) <p&(r 2 ). 


It is well-known that for each i, these problems admit a sequence of eigenvalues with 


Pi < P 2 < ■ ■ ■ < Pj < ■ ■ ■ 

Therefore, if we denote by (<£r)* the eigenfunction associated to (3j , we obtain that ip} = 
(+RYj{+p)i is an eigenfunction for the mixed eigenvalue problem (3) with eigenvalue 0). More¬ 
over, from previous expression of the eigenfunctions and the dimension of 'H 1 ((ri,r 2 ) x P ), we 
deduce that all the eigenfunctions are obtained by this process, i.e., we have proved 

Theorem 3.2. Let ip be an eigenfunction of ([3]) with associated eigenvalue (3. Then ip can be 
written as the product of two functions </?r : (ri,^) —» M and p>p : P -» M where ipp is an 
eigenfunction of Ap with associated eigenvalue (3, and v?r is a solution of the following problem: 


1 


1 


n— 1 


a 


d r ( a n 1 d r ip R ) + ^P+R = (G + 0) (p®., 


( 8 ) 


-^r(»t) = (J(n) + p) <£>rOt), 

Pr(t 2 ) = (J (r 2 ) + P) <PM.(r 2 ). 


The remainder of the section is devoted to obtain some additional information about the 
eigenvalues of this mixed problem on the warped case. As before, let us denote by {Pi} the 
ordered sequence of eigenvalues of the Laplacian Ap and by 0} the first eigenvalue of the 
Sturm-Liouville problem ([8]) with f3 = Pi. Then, the Rayleigh quotient applied to Problem ([8]) 
reads: 


-a n Vr^'r 


(9) p\ = min 

<PRe^ 1 (ri,r2) 


rr 2 


+ / a n 1 [O'r) 2 + (a 2 Pi - G)<£r] dr 


•r i 


rr 2 

/ a n "V B 

J n 


Tr dr 


With this characterization, we are able to prove the following two results: 
Lemma 3.3. For all i 6 N, 

P\ < P\ +1 - 


4 By notation, the dot will denote derivative of a real function. 
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Proof: For a fixed function the function on the right of Q is non-decreasing respect to 
Pi. Therefore, 


77 — 1 

-or L tpRipR 


P\ = min 

‘Pme'H . 1 (ri ,T2) 


r r 2 


+ / a n 1 [ Vr) 2 + (ot 2 Pi ~ G0</?r] dr 
r l Jr i 


a n Vr dr 


-a 11 Vm^'m 


r 2 r r 2 

+ / « n_1 [Vr) 2 + {or 2 Pi +1 - GVr] dr 

»T Jri 


nun 

<PR€'H 1 (ri,r2) 


' r\ 


rr 2 


a n l p\dr 


< 


= P\ +1 . 


’ r i 


□ 


Proposition 3.4. Consider the mixed eigenvalue problem ([3]) for a warped metric. If, for some 
*o £ N, it simultaneously satisfies J < 0 and a~ 2 Pi 0 > G for all r 6 (ri,r 2 ), then all the 
eigenvalues of the form /3* with i > ?'o are non-negative. If, in addition, at least one of the two 
inequalities is strict, then such eigenvalues are positive. 

Proof. Let {/3j°}jeN be the sequence of all the eigenvalues of the Sturm-Liouville problem ([8]) 
with P = Pi 0 . According to the first element of this sequence is characterized by 


-a n Vr^R 


P 




r 2 


rr 2 


+ 


= mm 

^RE‘H 1 (ri,r 2 ) 


r l 


u n 1 [(</%) 2 + (« 2 A 0 - G ) Vr] dr 


f r i 


rr 2 


a n Vr dr 


'r i 


Therefore, assuming that (fu reaches such a minimum and using the boundary conditions in 
0, yields 


Pi ( J « n Vr dr + a n 1 (VrV) + ^rV)) ) 

= —a n ~ l Vr( r 2 )j(r 2 ) + VrVVVi)) 

+ f ot n 1 [V’r) 2 + (a 2 Pi 0 — G) Vr] dr. 

J r\ 

As we can see, the right term of the above equality is non-negative under the hypothesis 
J < 0 and a~ 2 Pi 0 > G ; and it is positive if one of previous inequalities is strict. Therefore, p l ° 
is non-negative or positive accordingly. Finally, the result follows by recalling that 


Pi <K<% 


where we have used Lemma 3.3 for the first inequality. 


□ 


Remark 3.5. Previous result is optimal in the sense that, if we have that J = 0 and a~ 2 Pi 0 = G, 
then zero can be an eigenvalue of Q. In fact, if a = 1, the constant function is an eigenfunction 
of the zero eigenvalue. 









10 


S. C. GARCIA-MARTINEZ AND J. HERRERA 


As we will see in the forthcoming sections, Proposition |3.4| is particularly interesting in order 
to obtain results for both rigidity and bifurcation, whenever J < 0. If G < 0, it says that all the 
eigenvalues of the mixed problem are non-negative. If G > 0, it will restrict the possible zero 
eigenvalues to a finite family of Sturm-Liouville problems (see Section 5.2). 


4. Rigidity under metric variations 

4.1. Stating the variational problem. We will consider the variational approach first pre¬ 
sented by York [20 J , Gibbons and Hawking m, which has been further studied by Araujo in 
[4]. In these works, they consider a functional whose critical points are metrics with constant 
scalar curvature on the manifold M and where the boundary dM is composed by a minimal 
hypersurface Si and a CMC hypersurface S 2 . 

Let us begin by considering the so-called Gibbons-Hawking-York functional, GHY-functional 
for short (also known as the total scalar curvature plus total mean curvature functional): 


F : M 


Q 1 — >F(g)= R(g)dv(g) + 2 H(g)dcr(g), 


IM 


IdM 


where R(g) is the scalar curvature in M and H{g) is the mean curvature on dM. The first 
variation of such a functional reads as follows (see Section [hj Equation (35) in the Appendix for 
details) 


5F g {h) = - [ (Rij - \gijR{g)) h l3 dv{g) 

(10) Jm V 2 / 

- / (Hq/3 - H(g)g a p)h al3 d<j(g), 

J dM 

where R. % j , h lJ and IR^ are the components for the Ricci tensor, the (0, 2)-tensor h E T g M{= 
S k (M )) and the second fundamental form II respectively. From here, we deduce the following 
result 

Proposition 4.1. A Riemannian metric g is a critical point for the functional F if and only if 
the Riemannian manifold (M, g) is Ricci flat and dM is totally geodesic. 

Therefore, the GHY-functional defined in the whole set of Riemannian metrics does not 
determine the variational problem of our interest. In order to solve this, we will restrict the 
domain of the functional to a particular subset of Riemannian metrics. For the sake of clearness, 
such restriction will be performed in two steps. On the first one, we will consider Riemannian 
metrics satisfying the constraint C a j,(g) = 1, where: 


C a ,b(g) = a dv (g) + b da(g 2 ), 

Jm J e 2 

with a, b real numbers, a > 0, and g 2 denoting the metric g restricted to S 2 . 

Remark 4.2. 
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(i) It is clear that the above constraint was deduced from (jT]) , since in our problem only X 2 
must have constant mean curvature, and not all the boundary. In fact, as it is known, 
f M dv(g) determines the volume of M while da(g 2 ) determines the area of X 2 . 

(ii) Note that if a = 0 (resp. 6 = 0) then 6=1/ da(g 2 ) > 0 (resp. a = 1/ J M dv(g) > 0). 
So, it is not restrictive to assume that if a = 0 then 6=1; or, in other case, if 6 = 0 then 
a = 1. 

Let us denote by M. a ,b the space of metrics under the following constraint 

:={g <E M : C a)b (g) = 1}. 

In order to study the critical points of the GHY-functional restricted to M. a ,b, the Lagrange 
multipliers method leads us to the study of the following functional 

(11) F X (g) = F{g) - A (Ca tb (g) ~ 1), for some A 6 M, 
whose first variation takes the following form 

6 {F x ) g {h) = [r, 13 - (fl( g) + A °) g ^ h ii d v(g)+ 

+ (Ha /3 ~ H(g 1 )g a0 )h a,3 da(g 1 )+ 

+ J^ (i l a p ~ {2H(9 ^ + Ab) <m) h a ?da(g 2 ) 

where g k denotes the metric g restricted to Then, by a classical argument involving the first 
Bianchi identity, we obtain the following result 

Proposition 4.3. A Riemannian metric g on M n is critical for F restricted to the space M a ,b 
if and only if g is Einstein, Xi is totally geodesic, X 2 is totally umbilical with constant mean 
curvature and, if R(g) denotes the scalar curvature and H(g 2 ) the mean curvature in X 2 , then 

(12) (n-l)bR(g) = 2naH(g 2 ). 

In the second step, we restrict the domain of F to the space of conformal metrics of g which 
also lie in M a ,b-, that is, we will consider the subset of Riemannian metrics 

Conf 0j ft(p) := M atb fl Conf(y) 

where Conf(g) is defined in (11- 

Remark 4.4. Recall that our aim is to use the fiber bundle version of the Implicit Theorem, 
and so, a Banach space as Fl 1 will be enough for our purposes. However, for the study of 
bifurcation points, we have to deal with the fiber bundle version of the classical result given by 
Smoller and Wasserman for bifurcation ([23]. This result requires additional technical conditions 
as Fredholmness, Palais-Smale, among others, and so, the Sobolev space has not the required 
regularity. In such a case, it is necessary to restrict even more the conformal class to the Holder 
space C k,a with k > 0 and 0 < a < 1 (see [16] for detailed studies). 
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The set Conf a) &(g) is a smooth submanifold of Conf(g) because it is the set of regular points 
of the map C a ,b{g), being its tangent space T g (Confab (g)) identified with the set 

(13) H\ b { M ) = iyf G H l (M) : / > 0 and Y j M f dv (d) + ^ 2 ~ £ fdv(g 2 ) = °| • 


Thus, under this last restriction, the first variation of the functional F x becomes (see (36), (37)) 




IM 


R(g) - n 


( R(g ) + A a) 


fdv(g)+ 


- (n - 2)H(g 1 )fda(g 1 )+ 
J Ei 


'S 2 


A 6 


(n - 2)H(g 2 ) - (n - 1)— fda(g a ) 


with / G Therefore, we finally conclude that 

Proposition 4.5. 4 Riemannian metric g is a critical point of the functional F restricted to 
the space Confab if and only if the scalar curvature R(g) is constant on M, Si is a minimal 
hypersurface (H(g l ) = 0), and S 2 has constant mean curvature H(g 2 ). Under these assumptions, 
R(g) and H(g 2 ) are related by (12). 


Remark 4.6. For prescribed R(g) and H(g 2 ), the values of a and b are determined by the 


equations C a fi(g) = 1 and (12) (see also Remark 4.2). Moreover, if g is one of the critical points 
described in Proposition 4.5 it will be also critical for the functional F x defined in (11) with A 


solving one of the following equations (if a 0 7 ^ b, both equations define the same A assuming 

©): 


(14) 


b\ = H(g 2 ) 2 ( n ^ , and a\ = R{g) U 2 


n — 1 


n 


This last restriction lead us to a variational problem whose cr itica l points satisfy the desired 
properties. As we brought to remembrance at the end of Section [ 2 . 1 | the criteria for rigidity (as 
well as for bifurcation) make use of the second variation of the functional over critical points. 
For g a critical point of F X . with A depending on both, R(g ) and H(g 2 ), the quadratic form 
associated of such second variation takes the following form on 'Hf a b (M) (see (38), (39) and the 
Appendix for details) 


6 2 (F x ) g (fJ) = 


(n — 2 )(n — 1) 


(n — 2 )(n — 1) 


1M 


>M 


\grad(f)\ 2 g - 


/A s /- 


R(g) 

{n~ 1) 

R(g) 


f 2 


dv(g ) - 


IdM 


H(g) 
(n~ 1) 


fda{g) 


(n- 1) 


r 


dv(g) 


IdM 


-fdnf - 


H{g) 

(n~l) 


r 


da(g) 


where \grad(f)\ g denotes the squared norm of the gradient of /, A g is the Laplace-Beltrami 
operator and d n is the inward normal derivative (recall that, for the second equality, we have 
used the first Green identity). Then, we can recover the expression for the second variation 
and, even more, we can describe it in terms of Fredholm operators by using the following inner 
product on the space L 2 (M ) n L 2 {dM) 
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S(k,f)= / kf dv(g) + / kfda(g). 

JM JdM 

In particular, the second variation takes the following form 


6 2 (F x ) g (f,k) = 


(n — 2 )(n — 1) 


£(Jg(f),k) 


where J„ is a linear elliptic operator given by 


'Ig | M — Ag 


R(g) 

(n-iy 


jq | 8 M — d n 


H{g) 

(n- 1)' 


From the second Green identity follows that J g is a self-adjoint operator relative to the L 2 (M) n 
L 2 (dM)~ inner product. Namely 


£(Uf),k) = j M (kA g f - kfj *( 9 ) - J gM (»„/ + Ms) 

= L ( /A,k “ (Y^ fc/ ) dv(a) ~ L ( /fW + ^T) kf ) Ms) 

= £(f,J t (k)). 


Therefore, J g coincide with the Jacobi operator associated to 5 2 {F x ) g . Then, for our studies 
on rigidity, we have to study the spectrum of J g which lead us to the mixed eigenvalue problem 


(15) 


A q f - R ^ 9 \ f = Jlf in M, 

9J (n-\) J hJ 

-dnf ~ 7—= Pf on dM 
(n- 1) 


with / € H l a b {M). 

Remark 4.7. It is worth pointing out that in spite of the general case studied in Section |3j 
here the eigenfunctions should belong to This means that not all the eigenvalues of 

previous mixed problem are considered, but instead a subfamily of them. We will call admissible 


to an eigenvalue of (15) whose associated eigenfunction satisfies the integral condition in (13). 


4.2. General rigidity results under metric variations. The lower bound for eigenvalues 
obtained in Proposition 3.1 allows us to give a general result for rigidity associated to metric 
variations by giving a very simple and geometrical argument. 


Theorem 4.8. Let M be a compact manifold with boundary dM = Si U £ 2 . Consider 
a family of metrics which are critical points for the functional J~ x restricted to Conf^ for some 
a = a{ A) and b = b{ A); that is, metrics with constant scalar curvature in M and with zero mean 
curvature on Si and constant mean curvature on £2 . If R(g\ t ), J^((5a,) 2 ) ^ 0 for some A* £ I, 
being one of these inequalities strict, then the family {gAlAe/ locally rigid at A*. 
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Proof. Let us study the spectrum of J gx for A = A*. As we have pointed out in the previous 
section, the spectrum of such a linear map is related with the mixed eigenvalue problem detailed 
on (15) with g = g\ t . Now, Proposition 3.1 ensures that all the eigenvalues for the Jacobi 
operator are positive and g\ t is a non-degenerate critical point. Then, the result follows from 
the fiber bundle version of the Implicit Function Theorem (see m Appendix A]). □ 


5. Metric rigidity for foliations in warped product spaces 

Warped product spaces determines a foliation of the entire space by CMC-hypersurfaces 
known as slices (the natural leafs of the product). Such a foliation can be related with a family 
of metrics defined on a (fixed) compact subspace with the structure of a manifold with boundary. 
In particular, under the additional assumptions of constant scalar curvature on the space and 
the minimality of one of the leafs, our previous results on metric rigidity are applicable. As the 
construction of the metrics from the foliation is reversible, such a rigidity will be interpretable 
as metric rigidity for the natural foliation on warped spaces. 

Let us consider 

(16) (M n , g) = ((ri,r 2 ) x P n_1 , dr 2 + a(r) 2 g p ) 

an n-dimensional warped product space with constant scalar curvature where P is an (n — 1)- 
dimensional closed (compact without boundary) manifold. For each 7 E [ri,r 2 ), let us denote 
by S 7 = {7} x P a leaf of the natural foliation by CMC-hypersurfaces of the warped product. 
Along this section, we will assume that the metric is extensible to the slice Si := S ri and Si 
is a minimal hypersurface. Moreover, we will also assume that r 2 6 lU {00} (the case with 
{r 2 } x P minimal is analogous). 

For 7 E (? 7 ,r 2 ), consider the slab hL, = [ 7 - 1 , 7 ] X P a compact manifold whose boundary is 
composed by the leafs Si and S 7 of the foliation (see Figure [2|, being the (constant) mean 
curvature of the latter given by 

(17) ff(7) = -("-l)4v 

«(7) 

(recall that, as Si is minimal, previous equation implies that a(r 1 ) = 0 ). 

By a standard procedure, for a fixed ?’o E ( 77 , r 2 ) we can define a diffeomorphism T 7 : fi 7 —> 
n ro preserving the orientation of 8/dr (for 7 = ro, such a diffeomorphism is just the identity). 
Then, we define a family of metrics {(? 7 } 7g ( T . 1]r2 ) 011 ^ r„ given by the push forward of the metric 
g of fl 7 on Q ro , i.e., 

(18) g 1 = (T 7 )* 5 . 


Remark 5.1. As it is clear from the above construction, both (fi ro , g 7 ) and (fi 7 , g) are isometric 
under J/ 7 , and then, they share the same geometrical properties. So that, we can make the 
computations with the latter, which will be simpler in practical cases. Moreover, this approach 
shows the independence on ro of our results. 


By taking into account the properties described for the slab I) 7 and previous remark, we 
have that g 7 is a critical point for the functional being A : ( 77 , r 2 ) ->la function given 

implicitly on (14) (recall that a and b are determined by R(g^) and R((<? 7 ) 2 ) = H( 7 ), see 
Remark 


4.6). So, we are ready to apply Theorem 


4.8 


to the family {.9 7 } 7 e( ri ,r 2 ) defined on fi ro . 
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f 


P 


n —1 


d 




Figure 2. Construction of fi 7 ’s. 


However, it is interesting to present the result on this section, not as a property for the family 
of metrics, but for the family of compact sets {^ 7 } 7 e(ri,r 2 )- I n this sense, we will say that the 
family {H 7 } 7g ( ri r2 ) has a bifurcation point on 7 * if there exist sequences {/ n } n , {7 n }n with 
f n G 'H 1 (H 7n ), 7 n G (77, r 2 ) such that: 

( 1 ) lim 7 n = 7 * and lim g n = g Xt where g n := f n g. 

n—>00 n —>00 

( 2 ) For all n, the metric g n defined on fl 7ri determines constant scalar curvature R(g)', and 
makes Si = {n} x P a minimal hypersurface and S 7 = { 7 } x P a hypersurface with 
constant mean curvature 

(3) For all n, f n # 1. 

If 7 * is not a bifurcation point, we will say that the family {H 7 } 7 is locally metrically rigid at 

7*- 


As it is clear from construction, rigidity and bifurcation of the family {H 7 } 7 are equivalent 
to rigidity and bifurcation on the associated family of metrics {fl 7 } 7 . Then, Theorem 4.8 lead 
us to 


Corollary 5.2. Let ( M,g ) be a warped product space as in (16) with constant scalar curvature 
R(g ) and being Si := {77} x P a minimal hypersurface. Let us consider {^7} 7 e(n,r 2 ) a family 
of compact sets defined by H 7 = [77,7] X P. If R(g),H( 7*) < 0 for some 7* G (?7,r2), being 
at least one of these inequalities strict, then 7 * is a metrically rigid point for the family of sets 
{f2 7 } 7 . 


5.1. Some examples of rigid warped spaces. When ( M,g ) is a warped metric as described 
in (16), the scalar curvatures of both g and g p are related by 
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(19) 


R(g ) = -——- ( R(g P ) — ((n — 2)d 2 + 2da)) . 


As it is clear, whenever the warping function a and the scalar curvature of the manifold P 
R{g p ) are constant, the scalar curvature of the warped space R(g) is also constant. Moreo ver. 
all the slices in such a case are minimal, so H( 7 ) = 0. Therefore, if R(g p ) < 0 by Corollary 5.2 
we have that the family {n 7 } 7 is metrically rigid. 

However, we can obtain others non-trivial results. Let us assume that n = 3, R(g p ) is 
constant and R(g) is a non-positive constant. Then, by writing R(g) = —QE for some non¬ 
negative constant E, we can deduce 


( 20 ) 


2 k 

d(r ) 2 = R(g P ) -——- + Ea(r) 2 , whenever K to be a constant. 


Indeed, from 


a(r) 

we have that R(g p ) + 3 Ea 2 = d 2 + 2 da multiplying by d we obtain 
(. R(g P ) + 3 Ea 2 ) d = d 3 + 2do;a = (ad 2 ) 

and integrating we have 

R(g p )a + Ea 3 = ad 2 + 2K, 
where K is an integration constant, and ( |20[ ) follows. 

Next, taking the variable change 


( 21 ) 

the warped space becomes 


s = a(r) 


( 22 ) 


(s, s 2 ) X P, 


R(g p ) ~'^+Es 2 


: ds 2 + s 2 g p 


where s'> 0 is a zero for R(g p ) — 2 K/s + Es 2 and s 2 G M U { 00 }. As we can see from © and 
( 21 ), the mean curvature function IL(y) takes the form 


H( 7 ) = -2 


VR(g p ) - 2K/7 + 


7 


So, it is negative except at 7 = s where is zero. Thus, the metric space described in ( 22 ) 
satisfies that R(g),H( 7 ) < 0, being the inequality for the mean curvature strict for 7 > s'. In 
conclusion, Corollary |5.2| lead us to 


Theorem 5.3. Let (M 3 ,g) be a warped metric as in (16) with fiber of constant scalar curvature. 
If R(g) is a non-positive constant, a is monotonically increasing and d[r\) = 0, then the family 
of compact sets {H 7 } 7 is metrically rigid. 


There are several examples with R(g p ) constant. For instance, if we look for R{g p ) = —1, we 
can take P = M 2 /T where H 2 is the 2-dimensional hyperbolic space and T a co-compact group 
of isometries defined on it. For R{g p ) = 0, we can take P = T 2 the 2-dimensional flat torus. 
The case for R(g p ) = 1, where P is naturally the 2-sphere § 2 , is specially interesting due to its 
relation with general relativity. In fact, for E > 0, the space 
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(23) 


((s,Oo) X § 2 ,gK,E) , 


9k,e ■= 


1-^ + Es 2 

S 


i 2 , 2 S 2 

zds + s g 


represents the spatial slice of the Anti-de Sitter Schwarzschild spacetime. The constant K de¬ 
notes here the mass of the system, and so, it is considered positive. E is the so-called cosmological 
constant. In the particular case where E = 0, we obtain the Schwarzschild model. 

As we can see in {251 Chapter 6], the Schwarzschild models appear as a solution of the Einstein 
equations under some mild hypotheses. On the one hand, it is considered the case of vacuum 
Einstein equations, which lead us to the so-called Einstein metrics (metrics whose Ricci tensor 
is proportional to themselves). On the other hand, it is assumed that the spacetime is static. 

The minimal fibre Si = {s} x § 2 has a nice interpretation from the viewpoint of relativity. It 
determines the so-called marginally outer trapped surface (MOTS for short) which represents 
the boundary of the black hole contained in the model. That is, Si is a two-sided surface where 
the outgoing lightlike rays emanating from it collapse marginally (see [3] and the references 
therein for general background on MOTS). 

In particular, the previous results and the definition of metric rigidity allow us to deduce that 
small variations of the metric inside the conformal class, fixing the scalar curvature R(9k,e) and 
making Si minimal cannot have slices {7} x § 2 with constant mean curvature H( 7). In the 
context of relativity, the results imply that, once it is fixed the causal structure of the model 
(which is determined by the conformal class of the metric), the spatial fibers of Anti-de Sitter 
Schwarzschild spacetime are locally unique in the family of Einstein metrics preserving MOTS 
and making {7} x § 2 a CMC surface with prescribed mean curvature. 


5.2. Warped product manifolds with positive constant scalar curvature. Notice, in the 
previous rigidity results, we always require that the mean curvature and the scalar curvature are 
non-positive. However, if we look at Proposition |3.4[ the result gives some additional information 
in cases where the scalar curvature is positive. In this last section, we will show a simple method 
to obtain both rigidity and bifurcation results, when the mean curvature is negative but the scalar 
curvature is positive. 


Let us consider (M, g) a warped metric space as in (16) and assume that R(g) is a positive 
constant. We will also assume that the warping function a is non-decreasing (i.e. a > 0), being 
zero in 77 . In particular, and recalling the expression for iL(y) in ( |17[ ), the mean curvature of 
the slice {7} x P is nonpositive for 7 £ ( 77 , 77 ) and zero in 7 = 7 - 1 . Oru' aim is to study the rigid 
character of the family of compact sets {0 7 } 7 and, to do so, we have to study for each 7 the 
spectrum of the mixed eigenvalue problem ( |15[ ) with M = fi 7 . 

Let us recall two important facts. First, that the eigenfunctions / of previous mixed problem 
should belong to the space 'H\ b ( H 7 ), where a and b are determined by R(g) and H( 7 ) (Remark 
4.6). Secondly, and recalling Theorem 3.2, that any eigenfunction / of the problem can be 
written as the product of two functions fp and /®, where the former is an eigenfunction of the 
Laplacian on ( P,g p ) and the latter is an eigenfunction of the following Sturm-Liouville problem 
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(« n ] d r f M ) + ^ 


(24) 


= (R(g) + y) h, 


= yfa(ri), 


Here, /x denotes the eigenvalue associated to fp. T he e igenvalues of previous Sturm-Liouville 
problem will be denoted, in accordance with Section 3.1 as {p)( 7 )}. 

As / = /r/p G 77* 6 (fi 7 ), we deduce that: 


(25) 


( na 

\J 


a 


n —1 


(r 


r)dr+ ^ U 0 1 ^ /m( l)ot n X (t) 


'r 1 


fpd&p = 0 , 


where dcrp denotes the volume form on ( P,g p ). The second integral of previous expression is 
zero if and only if fp is an eigenfunction associated to fM with z 7 ^ 0. When i = 0, a restriction on 
the function /m appears. Therefore, any eigenvalue for (24) with n = m and i > 1 is admissible 


in the sense described in Remark 4.7 However, when fi = 0, we have to check additionally that 
the associated eigenfunction /m satisfies an integral condition derived from (25). 


Let us assume that a is bounded from above. Then, taking into account that linij^oo/LZj = 00 
and a is non-decreasing, we can find zq such that for z > zq, 


yiOt{r ) 2 > iMa(r 2 ) > R{g ) 


(here, we denote by a(r 2 ) = lim r _> r2 a(r)). So, as iL(y) < 0, Proposition 3.4 ensures that all the 
eigenvalues (p* ( 7 )} with z > zo are positive. Therefore, we only have to study when {p)( 7 )} 
is zero for 0 < z < zo- This study can be performed by using numerical approximations in the 
following way: For each z = 0,1,... ,zo and any 7 G {ri,r 2 ), consider the differential equation 
(24) with Jl = 0 and initial condition^] /r ( 7 ) = 1 and /m( 7 ) = H{^)- Then, two possibilities 


arise: 


(i) /R(ri) / 0, and then, /r does not satisfy the initial condition in (24) with ~p = 0. In 
particular, we deduce that zero is not an eigenvalue in this case. 

(ii) /R(ri) = 0, and so, zero is an eigenvalue for the mixed problem. In this case, and using 
the continuous dependence of the eigenvalues regarding 7 , we can make a local numerical 
analysis to determine if such an eigenvalue changes its sign around 7 . 

If for all 0 < z < zo we are in the first case, we have that 7 is a (local) rigid point. Otherwise, 
we have to analyze the cases where the zero eigenvalue appears in order to determine if the Morse 
index varies, and so, if 7 is a bifurcation point (recall Remark |4.4| about bifurcation results). 
In this second case, we have to take special attention when z = 0 in order to ensure that the 
eigenvalue is admissible. 


6. Appendix 

In this subsection we will include all the basic computations needed in order to compute the 
first and second variation of the functional J rX . Such computations have been already appeared 
elsewhere (see for instance HU SI), but we will include them here for the sake of completeness. 

^Observe that fmi'y) cannot be zero as, otherwise, /r( 7 ) = 0 and so fm is constantly zero. 
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Even if we are interested only in conformal variations, the computation of the first variation 
of the different elements conforming will be done for a more general family of variations. 
Let S k (M) the space of all symmetric (0,2) tensors of class C k {M) with k > 2. We consider a 
metric variation g : M k (M) x (—e, e) — > M defined by: 

g(h,t ) = g + th 

where M k (M) is the open cone of S k (M ) consisting of all Riemannian metrics on M such that 
for all g G M k (M) the tangent space T g M. k (M) is identified with the Banach space S k (M). By 
compactness, for |t| sufficiently small, g(h,t ) is in M k (M). 

Convention 6.1. Henceforth, all the elements associated to the metric g(h,t ) will be denoted 
as functions of t, assuming that the metric variation h is fixed from the beginning. So, elements 
as the metric itself, the scalar and mean curvature (among others) will be denoted by g(t),R(t ) 
and H (t) respectively. 

For the first variation of C a ^, just recall that 

(26) 5dv(t) = ^g lJ hijdv(g), 
and so 

(27) dC ab {t) = a /" 9 ZJh ij dv(g) + ^ [ g al3 h a/3 dcr(g 2 ). 

z Jm z Jt.2 

Now, we will focus on the first variation of the GHY-functional (recall the definition on (| 10[) ). 
From basic computations, we have that 


5F(t)= [ 6 (R(t)dv(t)) + [ 5 (H(t)da(t)) 

Jm JdM 

(28) = f (5 R(t))da(g) + R(g)(5 dv(t))+ 

Jm 

+ [ ( 6H(t))da(g) + H(g)(5da(t )) 

JdM 

and so, we need to compute the first variation of both the scalar curvature and the mean 
curvature. For the first one, let us consider a point pgM and a local normal coordinate system 
(xf) centred on p. Recall that in such coordinates, at the point p we have: 

9ij = Sij, d k gij = 0 and V* = d * 

for 1 < i, j, k < n, where 6^ represents the Kronecker delta, that is, Sij = 1 if i = j and Sij = 0 
if i j. In particular, the variation of the Christoffel symbols become^] 

(29) <Sr£-(i) = l -g kl (Vjhu + Vihjt - Vihij). 

^Observe that the variation of the Levi-Civita connection <5V is a tensor (in spite of what happens with V), 
and so, a special coordinate system can be considered for its computation. 
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From the definition of the scalar curvature, we have that: 

5R(t) = (Sg^WRij + g ii (SR ij (t)). 

In one hand, one easily obtain that 

d gij(t') = hij and 8g lJ (t) = —h lJ . 

On the other hand, the variation of the Ricci tensor has the following expression 

SRij{t)= 6Vi r^.(i) - <5VirJ. z (t) + <s (r?.r{ u - r“ r^) (t) 

= VidT^-VidT^t) 

— 2@ )V jhim “I - V'/V ihjm V/V mhij ViVjdl; m ) 

= V, (V,- V,/;') 


where we have used that 5 


pw pZ _ pu pZ 

1 ij L lu 1 Ij 1 iu 


(t) = 0, as the Christoffel symbols vanish for t = 0, 

the formulae (29) and the fact that g lm (S7ihj m — V m hji) = 0 for 1 < l,m < n. In particular, 
the variation of the scalar curvature takes the following form: 


(30) 


5R(t) = —h lJ Rij + Vi [Vjh' 3 - Vih\ 


Summarizing, the first two elements on the right of (28) bearing in mind (26) and (30), we 
obtain: 


(8 R(t))dv(g) + R(g){8 dv(t )) = 


IM 


im 


—h lJ Rij + -g v hij R(g) ) dv(g) 


+ / Vj (Vjti* - V*C) dv(g) 


IM 


or, by using the Divergence theorem, 


(31) 


8 / R(t)dv(t) = — 

Jm 


M 


/ dM 


h' 3 Rij ~ -g l3 h ij R(g)j dv{g) 
N l (V j hij - Vih[) da(g) 


where N denotes the inward normal vector to the hypersurface dM. 

Next, we need to compute the second term in (28). As the point p G dM, let us consider 
a different coordinate system. Take (x\,.. ., x n _i) a normal coordinate system associated to 
the boundary dM, endowed with the metric induced by g, and denote by 'yt Xl ,...,x n ) the geodesic 
starting at (xi,..., x n _i) with direction N. Observe that (xi,..., x n ) := ' y( X i,...,x n -i)( x n ) defines 
a coordinate system for M around p in which the metric g takes the following form: 

g = dx^ + g a pdx a dxP. 

Moreover, as (x a ) is a normal coordinate system for dM around p and N = d n , 


Kg = 
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Now, we are ready to compute the first variation for the mean curvature. Recall that, from 
definition, 

H(t)=g^(t)JL a g(t), 

and so, 


(32) 


5H(t) = (6g a P(m a p + g a P(8JL a p(t)). 


For the variation of the second fundamental form, define p : dM x (—e, e) — > ( TdM )-*- an unitary 
inward normal vector of the tangent space for the metric g(t). In particular, it is not restrictive 
to assume that, for p(0) = N = d n . Then, 


(33) 


8I a g(t) = 5gij(t)v*(t) T J a/3 

= h nn r^ + (Su(t)rr^ + sr^. 

Now, taking into account that Sr/(t)p*(t)p- , (i) = 1 and gi a (t )p*(f) = 0 for all t, we obtain that: 

(Su(t)) n = ~^h nn , 


which, together with (29), yields the following expression for (33): 

<5n apit) = \h nn T”p + \ (V a hp n + V ph an - V n h a p). 


Hence, (32) becomes 

5H(t) = -h^JL a p + \h nn g ^+ \g^ (V« hp n + V ph an - X7 n h a p) 


(34) 


-h^lap + l -h nn H(g) + V a h an - l -V n h\ 


a 

OL 5 


and so, using (34) and the fact that 5dcr(t) = (l/2)g a ^h a pda(g) we get 

S [ H(t)da(t) = f (~h a ^ a0 + \h nn H{g) + V a h an - \v n h% + ^H(g)g a ^h a p] da{g) 
J dM J dM \ Z 1 1 J 


J (-Ha /3 - \H{g)g a ^j h a0 + \h nn H(g) 

+ [ (v a h an ~l\7 n h«)da(g). 

J dM \ 1 J 


dcr(g) 


Then, using this identity, (31) and the fact that in our coordinates u = d n , we deduce that: 


5F= 5 R(t)dv(t) + 28 / H(t)da(t) 
Jm J dM 

h t3 Rij - \w 9ij R{g) ) dv(g) 


+ 


- (2H a/3 - H(g)g a p) h a ? + h nn H(g) da(g) + / ( V a h an ) da{g). 

J JdM 


J dM L 

Finally, taking into account that, 

V a h an = D a K + h a0 JL a p - h nn H(g), 
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where D denotes the induced connection on dM, and the fact that 



D a h%da(g) = 0 


from the Divergence theorem (recall that dM is closed), we finally obtain the following expression 
for the first variation of the GHY-functional 


(35) 


5F g (h) = - J ^ Rij - ^gijR{g)^j h l3 dv(g) 

~ [ (Ha/3 - H(g)g afj ) h af3 da(g). 

J dM 


Now, we will focus on the computation of the second variation. For this case, and for simplicity, 
only conformal variations will be considered (for the general case, see [3]). So, let us assume 
now that h = fg for a positive function / e T{ l (M). With this assumption, observe that (27) 
and (28) become 

(36) 

and 


«*</») = ?/■ /*(») +^ 


IM 


/S 2 


fda(g 2 ), 


(37) 


R(g) fdv(g) + (n- 2) j H(g) fda(g), 
z JM J dM 


2 JM 

respectively. Therefore, the second variation of the former is given by 


(38) 5 2 (C a , b ) g (fg) = 

where we have used that 


9 

n a 


fdv(g) + ^ ^ b I f 2 da(g 2 ), 

4 JM 4 Jy.2 


n 


Sdv(fg) = —fdv(g), 


71 — 1 

ddcr(fg) = —-— fda(g 2 ). 


For the second variation of the latter, observe that both (30) and (34) become 

1 71 — 1 

SR(fg) = -fR(g) + (n - 1)A gf, SH(fg) = — -fH(g) - —d n f, 


where A g = —div g (grad()) is the Laplace-Beltrami operator. 

Finally, the second variation of the functional F for conformal variations has the following 
form 


( 39 ) 


d 2 {F) g {fg) = j M f (^(sO (~9~) f + ( n ~ 1 ) a 9^’) dv (s) 

+( " “ 2) L 1 ( (m 1 ) h{3)/ ' (at) 8nS ) da(3) ' 
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